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I. Introduction 

The one-dimensional asymmetric simple exclusion process (ASEP) [131 El 
is one of the simplest interacting particle systems with a single conservation law 
(density of particles) and as such is a basic model in both probability theory and 
nonequilibrium statistical physics. Its importance is further enhanced by the fact that 
weakly asymmetric limits of ASEP distributions can be interpreted as distributions 
of the height function which solves the Kardar-Parisi-Zhang equation [121 El IT^ - 

Recall that in the asymmetric simple exclusion process particles are at sites of the 
lattice Each particle waits exponential time, then with probability p it moves one 
step to the right if the site is unoccupied, otherwise it stays put; and with probability 
q = 1 — p it moves one step to the left if the site is unoccupied, otherwise it stays put. 
Each particle does this independently of the other particles. For a finite number of 
particles this defines a Markov process; and for infinitely many particles, with further 
work [T3] this too defines a Markov process. 

In multispecies ASEP particles belong to different species, labelled 1,2,...,M. 
Particles of a higher species have priority over those of a lower specieso Thus, if 
a particle of species s tries to move to a neighboring site occupied by a particle of 
species s' it is blocked if s < s', but if s > s' the particles interchange positions. 
Second-class particles were introduced by Liggett [15] and subsequently developed 
and generalized by several authors [H HI [5l El [HI [HI [20] . 



^Many authors consider ASEP on the circle or the lattice [1, L] with open boundary conditions. 
^This is sometimes called the M + 1 species model, empty sites behaving as particles of another 
species. With our convention, a particle of species M is first-class, having priority over all others. 
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A configuration in ASEP with particles is the set of occupied sites 
X = {xi, . . . ,Xn}, {xi<---<xn)- 

Theorem 2.1 of [T8] (with proof corrected in [19]) was a formula for Py(A; t), the 
probability that the system is in configuration X at time t, given that the initial 
configuration was Y = {yi, . . . .Un}- It is a sum over the permuation group Sjq of 
multiple integrals. If p 7^ then 

^-(^' ^) = E 7^ Imo ne,) n (c-'^^e^^^'^*) d-e, (1) 

where Cr is a circle about zero in C with sufficiently small radius r, where 

and where Aa-{^i, . . . , C,n) is given explicitly by ([7]) below. 

In multispecies ASEP a configuration A* is a pair {X, vr) where X = {xi, . . . , Xjv} 
as before and tt is a function from [1, N] to [1, M]. If the system is in configuration 
X then the ith particle from the left is at Xi and belongs to species tTj. A special 
case is that of first- and second-class particles, a first-class particle having priority 
over a second-class particle. For example, if vr = (1 2 2 2) the left-most particle is 
second-class and the other three are first-class. 

The purpose of this paper is to establish for multispecies ASEP a formula anal- 
ogous to (lH) for Py{X; t), the probability that the system is in configuration X = 
{X, it) at time t, given that the initial configuration is 3^ = {Y, v). We show that 
there is an entirely analogous formula, 

but now the factors are not (except in special cases) given explicitly. They are 
determined by ([21]) and below. Note that they also depend on u. 

In the next section we present the proof of ([1]) in some detail. This is partly 
because for a correct proof one must refer to both [IB] and [11], but also because we 
shall show how to prove ([2]) by a modification of the proof of ([T]). 

There is another difference here. The factors must satisfy the family of iden- 
tities determined by (^T\i and (12 2 p below. For ordinary ASEP there is no issue about 
the existence of a solution since it can be written down explicitly. But for multi- 
species ASEP we must show that the identities define A'^ consistently. If, as in [1], 
the multispecies ASEP is formulated as a nested Bethe Ansatz problem, one is led 
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to show that the Yang-Baxter equations [6l [21] are satisfied. That the model has 
nontrivial solutions to the Yang-Baxter equations can be traced back to early work of 
Perk and Schultz on multistate vertex models [T6]. It was Alcaraz et al. [2] who rec- 
ognized that these Bethe Ansatz solvable multistate vertex models lead to integrable 
stochastic models. In our formulation below, these consistency conditions are stated 
in terms of representations of braid relations. 

In the last section we focus on the explicit determination of some A^, with par- 
ticular attention to the case of one second-class particle. 



II. ASEP with One Species 



A. Bethe Ansatz solution 



The probability Py{X; t) satisfies the differential equation (the master equation) 



du 
'di 



N 

\^u{xi - 1) (1 - 5{xi - Xi-i - 1)) + qu{xi + 1) (1 - 5{xi+i - Xi - 1)) 

i=l 



-pu{xi) (1 - d{xi+i - Xi - 1)) - qu{xi) (1 - 6{xi - Xi^x - 1)) 



(3) 



Here u is u[x\^ . . . ,xn', t), and in the zth summand above the jth variable is xj when 
j ^ i- We use the convention that any 6 term involving xq or xn+i is zero. 

The probability Py(X; t) is the solution of this equation that also satisfies the 
initial condition 

0) = 5y(X). (4) 

The particles interact through the exclusion constraint. If they did not interact 
then Py{X] t) would satisfy the differential equation 



du 



N 



= ^ \pu{xi - 1) +qu{xi + 1) - u{xi) 



(5) 



i=l 



If u satisfies this equation then it would also satisfy ([3]) if in addition the difference 
of the right sides were zero. This difference is 



N-l 



i=l 



^^(a^i, Xi+i - 1) + qu{xi + 1, Xi+i) - u{xi, x^+i) 6{xi+i - Xi - 1). 



Here we displayed entries i and z + 1 in u{x). For the first and last summands, we 
changed indices of summation from what they were in ([3]). 
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Thus u{x) satisfies (|3]) if it satisfies both ([5]) and, for i = 1, . . . , N — 1, the boundary 
conditions 

pu{xi,Xi) + qu{xi + + 1) - u{xi,Xi + 1) = 0. (6) 

For any nonzero complex numbers ^i, . . . ,^n, a solution of ([5]) is Yli {^i' e^*^^'^*). 
We may permute the ^j, take linear combinantions, and integrate. In this way we 
obtain a class of solutions 

where the functions Fcr{^) are arbitrary. We look for F^r such that the integrand 
satisfies the boundary conditions pointwise. This is the Bethe Ansatz. 

Substituting the left side of ([6]) into the part of the integrand that depends on Xi 
we get i^aii) Caii+i)^' times 

{P+ Q ^(7(i+l) - ia{i+l)) F^- 



Define Tjcr to be a with the entries a{i) and a{i + 1) interchanged. If we replace 
a by TjO" nothing else in the integrand changes, so a sufficient condition that the 
integrand is zero is that 

If we define 

Q(c c'\ - P + ~ ^ 

p + - ^ 

then the conditions become 

We can find the general solution of this system of equations for the F^j. An 
inversion in cr is a pair {i,j) with i > j and cr^^ii) < a^^{j). It is straightforward to 
check that one solution is 

^.(0= n (7) 

inversions 

All F^ are determined by F^,, where e is the identity permutation. Since this can be 
an arbitrary function (p{^), the general solution is F„{^) = Aa-{^) v^(0- 
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B. Satisfying the initial condition 



We choose so that the initial condition is satisfied by the a = e summand. 
Since 

where C is a circle about zero, we take = {2ni)^^ Yli^i^'^^ ^'^'^ domain of 
integration to be . 

When a 7^ e it matters which contours of integration we take because of the 
denominators in the A^j. If p 7^ then all the denominators will be nonzero on and 
inside Cr if r is small enough. It is such an r that we take in ([T]). 

Denote by I{a) the cr-summand in ([1]) with t = 0. To prove ([1]) we must show 
that 

5^/(a) = 0. 

This will be proved by induction on A^. For = 1 there is nothing to show. 
Assuming the result for A^ — 1, the sum over all permutations in iSAr\{e} such that 
(y{N) = N will be zero, by the induction hypothesis. (No ^-factor in A^r will involve 
^N, and there is an obvious correspondence between a G Sn-i and those a G Sn 
satisfying a{N) = N.) So it suffices to show that 

For each nonempty subset S of [1, A^ — 1] define 

Sj\f{B) = {(T G Sj\f : the inversions in a involving A^ are the (A^, i) with i G B}. 
We shall show that for each B we have 

Once we have this, ([H]) will follow since the set of a with cr{N) ^ N is the disjoint 
union of the various Siy{B). 

Start of the proof. The integrands in I (a) with a G Sn{B) may be written 

ieB i<N 
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In these integrals we make the substitution 



lU<n 

so that ?7 runs over a circle of radius . The integrand becomes 

i<N 

X n{^(^^'^^) ^ ^ > ^ > ^> < a-\£)}. (13) 

The reason we still have —1 in the exponents in f lT^ is that d^N = Y[i<N 
Lemma 1. When \B\ = 1 we have I{a) = for all a G Sn{B). 

Proof. There is a single i & B and (ITTl) is analytic inside the ^j-contour except for 
a simple pole at = 0. The exponent of in (fT2l) is positive since N > i, and so 
Vn > Hi, and o-^^(i) > a^^{N), and so Xo-i(j) > Xo-i(7v)- Therefore the integrand is 
analytic inside the ^j-contour, so the integral is zero. 

Lemma 2. When \B\ > 1 we have for all a G Sn{B), 

(id) 

where in the sum runs over all unordered pairs with i, j E B and i ^ j, where 
each Ib^ (ij) (cr) is a lower-order integral in which (fTT!) is replaced by a factor depending 
only on B and (the other factors remaining the same), and where = in the 
domain of integration. 

Proof. We may assume that q ^ 0. This case follows by a limiting argument. We 
are going to shrink some of the ^j-contours with i E B. Due to the defining property 
of r, the only poles we pass will come from the product f ITT]) . In fact, to avoid double 
poles later we take C,i G Cr^ with the all slightly different. 

Take j = max B and shrink the (^^-contour. The product (ITTl) has a simple pole at 
= (the j-factor has the pole and the ^-factors with i j are analytic there) and 
the power of in (fT2!) is positive as before, so the integrand is analytic at = 0. 
For each k E B with k ^ j we pass the pole at 
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coming from the A;-factor in (fTTj) . (Our assumption on the assures that there are 
no double poles.) For the residue we replace the A;-factor by 

P + QV n^^fc, TV - V Ue^N 

where in this and the j-factor we replace by the right side of ( 1141) . When i ^ j,k 
the 2-factor becomes 

and we replace in the numerator by the right side of ( fT4l) . 

We now shrink the (^^-contour. There is a pole of order 2 at .^a: = coming from 
( IT5|) and the j-factor in (fTTj) . Since k < j = maxB < N, we have Un — Uk > 2, so 
the exponent of C,k in ( fT2l) is at least 2. Therefore the integrand is analytic at C,k = 0. 
The factor ( JT5l) has no other poles inside Cr^. . An z-factor with i j,k will have a 
pole at ^k = C,i if fi < rk- There is also the pole at 

Kk — 7 

-p 

coming from the j-factor. This relation and f lT^ imply = C,k- 

Thus when we shrink the (^j-contour and the ^^-contours with k ^ j we obtain 
(A^ — 2)-dimensional integrals in each of which two of the .^-variables corresponding 
to indices in B are equal. This proves the lemma. 

Lemma 3. In the notation of Lemma 2, for each there is a partition of Sn{B) 
into pairs a, a' such that lB,{i,j){cr) + lB,{i,j){,0'') = for each pair. 

Proof. We pair a and a' if they agree except for the positions of i and j, which are 
interchanged. The factor f|T2|) is clearly the same for both when = ^j, and we shall 
show that the a- and cr'-factors in ([T3|l are negatives of each other when = ^j. 

Assume for definiteness that 

i<j and a'^{i) < a'^{j). (16) 

(Otherwise we reverse the roles of a and a'.) Then the factor S{^j, ^i) does not appear 
for a in f|T3l) but it does appear for a' . This factor equals —1 when = ^j. 

To complete the proof it is enough to show that for any k ^ i,j the product of 
S'-factors involving k and either i or j is the same for a and a' when = ^j. If (7~^{k) 
is outside the interval {a~^{i), cr~^{j)) the S'-factors in question are the same for a 
and a', so we assume a~^{k) is inside the interval. There are three cases, with the 
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results displayed in the table below. The first column gives the position of k relative 
to i and j, the second column gives the product of S'-factors involving k and either i 
or j for a, and the third column gives the corresponding product for a'. 



k<i S{^i,^k) S{^j,^k) 

k>3 S{^k.Xr) 

In all cases but the second the iS-factors are exactly the same for a and a' when 
= ^j. For the second we use S{^, ^k) S{^k, = 1- 

Clearly iQ follows from Lemmas 1-3, and this completes the proof of ([1]). 



III. ASEP with Multiple Species 



A. Bethe Ansatz solution 



Observe that an interchange of particles at positions Xi and Xj+i has the same 
effect as leaving the particles as they were but interchanging vTj and vFj+i. Thus X 
remains the same but n is replaced by Tin. (This is the same Tj as before, but applied 
to TT rather than a.) 

Write u'"{X; t) for Py{X; t). The master equation for w'^ differs from equation ([3]) 
for u since particles are not blocked as much, so there are other terms on the right 
side. 

We define 

{0 if TTj = TTj+i 
P if TTj < TTj+i 
q if TTi > TTj+i, 

and define /3i(7r) as above but with p and q interchanged. (Thus /3i(7r) = ai(Tj7r).) 
We compute that what must be added to the right side of ([3]) is 

N-l 

^ [ai(7r) u^'^'ixi, Xi+i) - /3i{n) u'^{xi, Xi+i) S{xi+i - - 1). 

i=l 

(As before, we displayed entries i and i + 1.) The terms here may be incorporated 
into the boundary conditions, as the 5-terms in ([3]) were. We conclude that if u'^{X; t) 
satisfies the equation 



dt . ^ 

1=1 



N 

J2 [pu^Xi - 1) + g^i"(x, + 1) - u"(x,)] , (17) 
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the boundary conditions for z = 1 , . . . , — 1 

pu^{xi, Xi) + qu^{xi + 1, Xj + 1) — u'^{xi, Xi + 1) 

- ai(7r) Xi + l) + A(7r) ^"(a;^, Xi + l) = 0, (18) 

and the initial condition 

u^{X-t) = 5Y{X)5M, 

then P3;(A'; t) = M^(X;f). 

We assume a solution of the form 

where = 61,(71) so that the initial condition is satisfied by the a = e summand. 

Substituting the left side of f lTS]) into the part of the integrand that depends on 
Xi we now get (^^(i) ^^(1+1)^' times 

[f{b,a) + (3,{n)b]Al-a,{7r)bAl^\ 

where 

f{^,e)=P + q^e-C, a = ^^), b = Caii+i). (19) 

If we replace a by Tjcr nothing else in the integrand changes, so a sufficient condition 
that the integral be zero is that 

[f{b, a) + A(7r) b] a: + [/(a, b) + A(7r) a] A^, - a.(7r) 6 A^^^ - a.(7r) a = 0. 

If this is to hold for all tt it must hold for tt replaced by Tjvr, so 

[f{b, a) + «,(7r) 6] Al'- + [/(a, 6) + a,,(7r) a] - A(7r) 6 - A(^) a^J,. = 0. 

If TTj+i = TTj then we obtain simply 

f{b,a)Al + f{a,b)Al^ = 0. (20) 

If TTj+i 7^ TTj we use /3j(7r) = 1 — Q;j(7r) to help with the computation and find by 
eliminating A^^ from the two equations above that 



b — a 



f{a,b) L f{a,b) 



a — b 



a: 



(Here (fT9l) was also used.) If we replace the second ai^ir) by 1 — /3i(7r) we get 

f{a,b) f{a,b) 
9 



which now agrees with (120|) when 7r.j+i = tTj. 

The relation is nicer for the quantity defined by 

Al = KA„, (21) 

where A„ is defined by ([7]) as before. The condition = 6^{7r) becomes h'^ = 5y(7r). 

If we use a — b = f{b, a) — f{a, b), divide by At^u = S{b, a) A^j, and recall (fT9ll . 
our formula becomes 

= K+{1 + [«.(vr) h^'"^ - A(7r) K]. (22) 

We have to show that these formulas together with h'^ = 6^{7t) define h'^ consis- 
tently. This means that if we obtain a from e by a sequence of operations Tj on 
then what we obtain from this same sequence acting on the h'^ is independent of the 
particular sequence chosen. This follows from the relations 

(i) Ti Ti = /; (ii) Tj Tj = Tj Ti when \i - j\ > 1; (iii) Ti Ti+i Ti = T^+i Ti T^+i. 

We should be more precise as to what exactly these Ti are. Let "Ho be the set of 
all functions 

h : n ^ function of ^, 
and n = SN X Uq. Define : U ^ Uq hj 

T^^a, h) = h+{l + 5(e.(o, e.(.+i))) W^ ■ {h o T,) - A ■ /^]. (23) 
Then define Tj : "H — "H by 

T,{a,h) = {T,aXi(^,h)). (24) 

(The same sympol Ti is used to denote these operators and operators on Sn- The 
context should make it clear which is meant.) It is the Tj defined by ( 124|) that 
satisfy relations (i)-(iii). Once we have these relations it follows from the Coxeter 
presentation of Sn [TT| Sec 1.9] that there is a homomorphism from Sn to the group of 
mappings from "H to itself, such that each Tj acting on iS^v goes to the corresponding 
Tj acting on "H. Because this is a homomorphism, if two products of transpositions 
in (Sat are equal then so are the corresponding products of the Tj on "H, which is what 
was claimed. 

Relations (i)-(iii) are the braid relations and are the consistency equations referred 
to in Section I. Relations (i) and (ii) can be checked by hand. Since (iii) involves only 
three consecutive indices, it is enough to check the case = 3. This was verified by 
a computer computation. 



10 



B. Satisfying the initial condition 

We have to show that 

E ^"(^)=0' (25) 

fTG5iv(B) 

where I^{cr) is the a-summand in ([2]) with t = 0. 

In Section II. B we made a variable change, shrank contours, and found that all 
I{a) with 0" e Sn{B) are sums of lower- dimensional integrals lB,(ij){cr), one for each 
unordered pair with i,j G B, such that = C,j in the domain of integration. 

To extend this to the we have to know first that no new poles arise from 

the factors h'^. These would come from S-factors involving C,n, because it is only 
these that might introduce poles after the variable change ffTOj) . Using fl22|) to see 
which such factors could arise in the expressions for h'^ we start with a = e and do 
the following to get to our a G Sn{B): first, bring to the front and rearrange the 
i ^ BU{N}. This gives some factors 1 + S'(^j, ^j) with i,j ^ N. (Which exact factors 
occur depends on vr as well as a.) Then move through the i E B. Since at each 
step was to the right of where it moves to, we get some factors 1 + S{^i, S,n) with 
i E B. Then rearange the elements of B to reach a, which introduces some factors 
1 + S{^i, with i, j G B. Thus the only factors involving are some 1 + S{^i, ^n) 
with i G Bu If we multiply this factor by the S{^N,^i) from the product in for 
A„ we get S{^iy,^i) + 1. Thus, no new poles arise from h'^ with a G Sn{B). 

We can now proceed as in Section II. B, making the variable change and shrinking 
contours. 

We say that a, a' are (z,j)-paired if they agree except for the positions of i and 
j, which are interchanged. We saw in the proof of Lemma 3 that if a and a' are 
(i,j)-paired then the integrands in lB,{ij){cr) and lB,{i,j){cr') are negatives of each 
other. Therefore to prove fl25l) it suffices to show that if a and a' are (z, j)-paired 
then h'^ = h'^, when = ^j. We show this by induction on \a^^{i) — a^^{j)\, and we 
may assume cr^^{j) > cr^^(i). 

If cr^^i^j) = (y^^{i) + 1 we replace i by cr~^{i) in fl22|) . Since S{^i,C,j) = —1 when 
6 = the statement holds then. 

Suppoes m > 1, that the statement holds when cr~^(j) = <7~^{i) +m — 1, and that 
(T~^{j) = cr'^ii) + rn. For convenience of notation we assume that 

a = (1, 2, m, m + 1 . . .), a' = (m + 1, 2, m, 1 . . .). 

•^For example, suppose cr = (3 2 5 4 1), for which B ~ {4, 1}. Then the steps might be 
(1 2 3 4 5) ^ (1 3 2 4 5) (3 1 2 4 5) ^ (3 2 1 4 5) (3 2 1 5 4) (3 2 5 1 4) ^ (3 2 5 4 1). 
The only S'-factors involving ^5 come from steps four and five, and are 5(^4,^5) and 5'(^i,^5). 
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(Thus a and a' are (1, m + l)-paired. Neither the actual labels nor the exact positions 
for the indices is relevant; only their relative positions is.) The dots represent other 
entries that are equal for a and a' . 

We want to show that = hi, when = ^m+i- We have 

TiO" = (2, 1, ... m, m + 1, . . .), Tia' = (2, m + 1, ... m, 1, . . .). 

Observe that Tia and Tia' are also (1, m + l)-paired but (ricr)^^(l) = 2 while 
{Tia)~^{m + 1) = m + 1, so the induction hypothesis holds. Thus k!^^^ = hl^^^i for all 
TT when ^1 = 

Applying (!22l) with 2 = 1 and a replaced by Tia and by Ticj' gives the relations 

K = K^^, + (1 + hl\l - A(7r) h-^J, 

K, = h\^, + (1 + 5(6, U^i)) [ai(7r) - /3i(7r) Z.^^,,]. 
When 6 = ^m+i the right sides are equal and therefore so are the left sides. 
This completes the proof. 

C. Formulas for 

1. A general formula 

For cr G (Sat we define G T^o to be the function tt — )■ /i^. With this notation 
formula (122|) may be written hT^a = ^i°(Ci ^o-), where T° is given by (!23|) . From this 

we get 

/iT.T,. = 7;°(T,a,/i^^,) = 7;0(T,a,T,°(a,/i.)) = T^T,{a,K). 
And in general, 

hr T ■■■T rr = T,„ ■ ■ ■ T,- (cT, 

If we first set a = e above and then choose the various Tj such that 

TnTj^---Tj^e = a, (26) 

we get, since /ig = 6^(7^), 

K = T^J,,...Tae,Su). (27) 
For p ^ Sn, and a and /3 functions of vr, we define 

[/,(p,/3) = l-[l + S(ep(.),ep(m) )]■/?, 
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Vi(p, «) = [! + S'(^p(i), ^p{t+i))] ■ a. 

These Ut and Vi come from the coefficients of h and h o Tj, respectively, in f l23|) . and 
products of them are the coefficients of the 5-summands which result from applying 
the various Tj in fl27|) . When we apply (1251) consecutively to flTTl) we get a sum of 
products of the form Wi W2 ■ ■ ■ Wm, where each W is a. U or V. 

As an example, we find that T° Ti{e, Si,) is a linear combination of Su, Si,ot,, S^oTj, 
and 5yoTioTj- The coefficient of 5„ is Uj(Ti e, f3j) Ui{e, the coefficient of S^oTi is 
t/j(rj e, f3j) Vi{ e, ai), the coefficient of 6uoTj is V^(7i e, aj) e, /^j o Tj), and the coef- 
ficient of S^oTioTj is ^^(21 e, ttj) Vi( e, o Tj). These are the four W-products. 

Passing to the general case, we observe that in every factor W the original per- 
mutation e has been composed with all earlier (rightward) Tj in ( 127|) . Every factor 
Vi comes from the corresponding hoTi summand in ( 123|) . This changes the a or /3 in 
each earlier factor W by composing it with these Tj. Otherwise said, each a or /3 in a 
factor is a composition with the later Tj. The VF-product is determined by these 
(ordered) i, and we denote the sequence of them by X = The resulting 

5-term is Suot,^o-ot,„- 

In general we get a linear combination of 5-terms of the form SuoTi-^o---oTi„, where 
X = {ii, ■ ■ ■ , in} C Xo- and Z„ = {ji, . . . ,jm} is the sequence in (l26ll . Each 4 is a 
j£ while the other je appear between some consecutive ik^i and ik. (For this we set 
^0 = 0, in+i = 00.) From the above remarks we see that the coefficient of (5i/ot,^o - ot,„ 
is 

m 

Wt = 1[Wxj (28) 

£=l 

where 

f ■ ■ ■ Tj^ e, a,, o ri,_^ o . . . o r,J if je = ik, 

Wx,t= I (29) 
I Uj,{Tj^_, ■ ■ ■ Tj^ e, /3j^ o Ti,_^ o ■ ■ ■ o J if G (zfc_i, 4). 

This gives 

/i:= 5^ Wi(7r)W,o...oT,„(7r), (30) 

Ida 

where X = {zi, . . . , 

The sequence 1^ in f l26|) is not unique, so we do not yet have an explicit formula 
for h'^. To find a particular sequence let us see how to get from cr to e by a sequence 
of Tj. We may do this by first bringing 1 to slot 1 by a sequence of transpositions, 
then bringing 2 to slot 2 by a sequence of transpositions, etc.. We first bring 1 to slot 
1 by the product 

Ti T2 • • ■ Tct-i(i)_i. 
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Suppose we have brought 1, . . . , — 1 to their slots. Then k itself has been moved to 
the right by 

L{k) := ^{j : j < k, > a^^{k)} = the number of inversions of the form {k,j), 

because bringing each such j to its slot has moved k one slot to the right. We then 
bring i to its slot by the product 

Tk - ■ ■ To.-i(fc)+t(fc)-i. 

Thus 

N-l 

Y\_^i Ta-^k)+L{k)-i cr = e, 

k=l 

the factor with the larger k being to the left. This gives 

N-l 

^ = Y\. '^<^-'^{k)+i-(k)-i ■ ■ - TkC, (31) 
k=l 

the factor with the larger k being to the right. With this representation we have 

N-l 

T.= [j {ik,---,k), (32) 

k=l 

where ik = cr^^{k) + L{k) — 1, the interval with the larger k being to the right. 
Together with ( 130|) this gives an explicit, albeit complicated, formula for /i^. 

2. One second-class particle 

Because all the tt now will have 1 in a single position and 2 in the others, most of 
the a and /? terms in ( l29l) will be zero when applied to n. So most of the Wi^ will be 
(when a Vj) or 1 (when a Uj). It is possible systematically to determine which Wx 
will be nonzero when X„ is given, and that makes it practical write down formulas 
for the /i^ in specific cases. We shall not describe the procedure but state two results 
that use ( 132|) . which comes from representation (13T]) . 

For z^~^(l) = 1 (the second-class particle initially in position 1) and 7r^^(l) = j 
(the second-class particle ending in position j) h'^ = when cr~^(l) < j and 

K={p-q S{Cu e.o))) q (1 + 5(6, ?.0-i))) ■■■q{l + >S(6, ?.(i))) , 
when o"^^(l) > j. (When (t"^(1) = j the factor on the left equals 1.) 
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For p ^(1) = 2 and vr ^(1) = j the formula is more complicated. When a ^(1) > j 
and CT~^(2) + i{2) > j it is 

j-i i-i 

i=l k=l 

i-1 

k=i+l 

fc=i 

We should point out that in any given case ( 131]) may not be the best representation 
for computation. For example, take z^~^(l) = 7r~^(l) = 3 and cr = (4 3 2 1). Using 
the representation 

cr = Tg Ta Ti Tg Ta Tg e (33) 

from ( 13 ip . there are five nonzero summands in ( 15U]1 . Therefore is given as a sum 
of five products. But if instead we use the representation 

a = TiT2TiTsT2Tie, (34) 

there are only two nonzero summands and we get the relatively simple formula 

K = q{l + 5(6, 6)) {q-p S{^2, 6)) q (1 + ^(6, 6)) 

+ (9-^^(6,6)) (p-g^(6,?4)) (g-P^(6,6)). 

For another example take a = (4 3 2 1) as before but z^~^(l) = 3, 7r^^(l) = 4. 
Using ( 15^ we again get five nonzero summands in ( 15U]1 . But using (IMj) we get only 
one, and the much simpler answer 

K = q{l + S{^^,^,)) (g-p5(6,6)). 

So although the representation ( 13T|) leads to a general procedure for computing 
h'^, in any particular case there may well be a better representation for computation. 
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